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ABSTRACT 




A study is made of the circumstances under which there will 
be incipient buckling in an infinite plate containing a circular 
hole and subjected to uniform uniaxial tension at a great distance 
from the hole. This buckling is a consequence of the fact that 
near the hole and in line with the direction of the applied tensile 
loading, there are areas in which both principal membrane stresses 
are compressive. Values of critical ten&Ll^ stress are given, 
together with expressions for the corresponding lateral deflection 
patterns. The results are in the form 

S = KE (t/a) 2 
cr 

where E is Young's modulus of elasticity, t is plate thickness, a 
is hole radius, and K is a buckling coefficient having the following 
values for the first four modes: 7.65, 7.80, 11.7, 12.3, A 

comparison is made with some experimental determinations of Danis 
for plates of finite width. 
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NOMENCLATURE 



a 


Hole radius 


A--- 

NM 


A set of unknown coefficients in a series 
describing plate deflection. 


A 


An unknown eigenvector 


A 


An integer index used in Appendix C 


a(J,L) 


A trigonometric product used as an intermediate 
function in Tables 5 and 6 


b 


The half-width of an infinite strip in convenient 
units 


b(J,L) 


A trigonometric product used as an intermediate 
function in Tables 5 and 6 


B 


A coefficient matrix in Appendix D 


B 


An integer index in Appendix C 


C 


An integer index in Appendix C 


d 


Hole diameter 


D 


An integer index in Appendix C 


D 


3 2 

Flexural rigidity of the plate. = Eh / 12 (1 - a) ) 


E 


Young's Modulus 


E^ ,E 2 ,etc . 


Terms in the membrane energy integral 


f 


Functions of I,J,k, and L used in Table 5 


g 


Functions of and L used in Table 6 


, G 2 , e tc . 


Terms in the bending energy integral 


h 


Plate thickness in units of hole radii 


I 


Identity matrix 


I 


A dummy integer index for the index N 


I' 


A single subscript representing the double subscript IJ 


I 


The maximum value of I* 


J 


A dummy integer index for the index M 
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NOMENCLATURE (continued) 



J' 


A single subscript representing the double subscript KL 


J 


The maximum value of J 1 


K 


A dummy integer index for the index N 


K 


Buckling coefficient based on hole radius 


L 


A dummy integer index for the index M 


M 


An integer index 


M 

r 


Bending moment per unit length 


M _ 
rt 


Twisting moment per unit length 


N 


An integer index 


N ,N 
x y 


Membrane forces per unit length, =(Th, etc. 

X 


P IJKL’ P I'J’ 


A term in the quadruple series that expresses the 
bending energy integral 


n 


A matrix of the eigenvectors of I 


q ijkl’ Q i'J' 


A term in the quadruple series that expresses the 
membrane energy integral 


r 


Distance from the hole center in units of hole radii; 
also a coordinate of the cylindrical system 


r ijkl’ R I'J' 


A term in the quadruple series that expresses the 
membrane energy integral 


R. 


The matrix of elements R 

IJKL 


S 1’ S 2’ S 3 


Functions of N, M and -b in Appendix A. 


S 


Externally applied stress 


s 

cr 


The value of S at which incipient buckling occurs 


t l ,t 2 ,t 3 


Functions of N, M and ^ in Appendix A. 


t 


Plate thickness 


T 


The matrix of elements T 

IJKL 


T T 

IJKL’ I f J' 


A term in the quadruple series that expresses the 
bending energy integral 


u 


Component of displacement in the x direction 
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NOMENCLATURE (continued) 



U l- U 2 



Membrane energy in the plate 



U 



3 



Bending energy in the plate 



v 



Component of displacement in the y direction 



W 



Plate deflection in units of hole radii 



W ,W 
r rr 



etc . 



x* 



x 



X NM 



Derivatives of plate deflection with respect to the 
subscripts 

A compacted notation of in Table 5 with I and J 

replacing N and M respectively 

A compacted notation for in Table 5 with K and L 

replacing N and M respectively; also a coordinate in 
the Cartesian system 

An unknown coefficient which is used to satisfy the 
boundary conditions for each term in N or M 



X 

y* 

y 



y nm 

cx 



An unknown eigenvector in Appendix D 

A compacted notation for in Table 5 with I and J 

replacing N and M respectively 

A compacted notation for Y in Table 5 with K and L 
replacing N and M respectively; also a coordinate of 
the Cartesian system; also a scalar parameter in an 
eigenvalue problem in Appendix D 

An unknown coefficient which is used to satisfy the 
boundary conditions for each term in N or M 

The number of terms taken in N in the series for W 



P 

x 

6 

t) 



The number of terms taken in M in the series for W 
A scalar parameter in eigenvalue problems 
Normal stress 
Shear stress 

An angular coordinate in the Cylindrical system 
Poisson* s ratio (*0 = .3 is used throughout) 



OTHER NOTATION 

A double underline indicates a square matrix. (A) 
A single underline indicates a column vector. (A) 
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1. INTRODUCTION 



The problem of determining the stress situation near a hole 
in a plate under tension (see Figure 1) has a well known solution 
which is of great practical importance. 

The original work on this problem was done by Kirsch (19), 
who solved the problem of the stress distribution in an infinite 
plate pierced with a round hole and subject to uniaxial tension 

i 

at a great distance from the hole. Howland (17) treated a similar 
problem with an infinite strip of finite width. Savin (5) and 
Hogan (15) have compiled solutions produced by other authors for 
variations in hole geometry, as square, triangular, etc,, for both 
finite and infinite plates. Lur'e (3) has solved the problem with 
a hole in a shell of uniform curvature. 

An important result obtained from these solutions is that the 
stresses near the hole become much greater than those imposed upon 
the plate at a distance from the hole. This has given rise to the 
wide publication of many stress concentration factors which are 
important in design work and which relate the maximum stress to the 
nominal stress at a distance from the hole. 

Another feature of the stress situation, which is not so widely 

recognized, is that certain areas of the plate are under uniaxial or 

k 

biaxial compression, although tension is imposed externally. These 
areas of compression, given the proper conditions of plate thickness 



* These terms will be employed herein to mean that one or both 
principal stresses are compressive. Specifically, "uniaxial compression 11 
means that one principal stress is compressive and one is tensile; 
"biaxial compression” means that both are compressive. 
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and external tensile load, will buckle, possibly causing failure of 
the plate. Houghton and Rothwell (16) have shown that in plates 
with holes having a very small gross curvature the stress concen- 
tration factor can increase appreciably. It is reasonable to assume 
then, that even though local distortions of a plate may be acceptable, 
such distortions could increase the stress concentration beyond that 
expected for a flat plate. Failure of the plate by yielding could 
then occur. 

There is no known investigation of this buckling problem, other 
than an experimental investigation by Danis (14), to determine the 
critical stress for the buckling of a finite plate with a round hole. 

This paper presents a theoretical investigation of the incipient 
buckling problem for an infinite plate containing a round hole. This 
problem is of interest in itself and for the practical reasons given 
in the preceding discussion. Results are presented on page 33. 
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2. PROBLEM DESCRIPTION 



Consider a plate infinitely extended in two directions and 
pierced with a hole, as shown in Figure 1. 

s 




s 

Figure 1 



x and y are Cartesian coordinates, r and 0, cylindrical coordinates. 
The z axis, not shown, is perpendicular to the paper at the center 
of the hole, and extends in the positive direction toward the reader. 
Consider all quantities to be in units of hole radius, a, which 
essentially reduces them to dimensionless quantities, with hole radius 
equal to one. 

The plate is placed in uniform tension, S psi, at infinity in 
the x direction. The solution of the stress distribution problem for 
this plate, as originally solved by Kirsch (19), and clearly presented 
by Timoshenko (7) is: 
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( 1 ) 







Y = — 
2 




( 2 ) 

( 3 ) 



This set of equations is henceforth referred to as the Kirsch 
solution. 

It is observed that the maximum tensile stress is 3 times the 
stress applied at infinity, and occurs at r = 1, and 9 = + 7T/2. 

The maximum compressive stress, equal in magnitude to the applied 
tensile stress, occurs at r = 1, and 9 = + 7f. From this it can 
be seen that if buckling were to occur in the plate, it would be 
most apparent near the hole along the x axis. 

A better understanding of the buckling problem is obtained by 
considering the minimum principal stresses (most negative, largest 
compressive). From the Kirsch solution, the minimum principal 
stresses are found and shown in Figures 2 and 3. 

Figure 2 shows one quarter of the plate, with the hole centered 
at the top left. The stress distribution is symmetrical with respect 
to the x and y axes. This picture was prepared with consideration 
only to whether the stress situation at a particular point is biaxial 
tensile, biaxial compressive, or uniaxial compressive. 

Figure 3 shows the same section of plate as Figure 2. This 
presentation maps the relative values of the algebraic minimum 
principal stress at each point. 
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FIG 3 - MAP OF THE MINIMUM PRINCIPAL STRESS 
AROUND A HOLE 
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It can be seen from these figures that there is an area having 
biaxial compression. Accordingly, buckling is expected to take place 



for some sufficiently great value of applied load. 

Numerous methods have been developed for solutions of plate 
buckling problems. In general, these methods involve either the 
satisfaction of the equilibrium equations or minimizing the strain 
energy of the plate (Energy Methods). A compilation of many solutions 
to plate buckling problems was made by Gerard and Becker (11). This 
is an excellent reference for both method and practical solutions. 

One of the first energy methods was developed by Ritz (20), and 
independently, by Lord Rayleigh (4). This method is presented in a 
rigorous manner by Sokolnikoff (6). Budiansky and Hu (12) give a 
detailed example of the application of this method to a plate buckling 
pr oblem. 

The Ray leigh-Ritz method involves determining the total strain 

energy in the plate as a function of an assumed buckled configuration. 

Equilibrium, and the stable buckled deflection, W(x,y), exists when 

this strain energy, , is minimized with respect to W(x,y). 

To do this, first, a family of functions,^ , are assumed, which 

are dependent upon parameters A , such that the buckled deflection, 

K 

W(x,y), can be represented by: 



The functions are substituted into the strain energy expression 
U so that 



W(*>3) = §(-x-,y,A„Ax, •' ■) 

where J is such that the boundary conditions on the plate are 



(4) 



satisfied for any values of the A. . 

l 




(5) 



where U is the integrated total energy in the plate. 
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In order to minimize the energy function, the following set 



of equations is formed 



« 

— - — o 

3 A, 

dU t 

d A z ~ ° 



( 6 ) 



If this set of equations has solutions A^,A , the function 



W('x ji/ )= §(x J y ) A h A X) t - • • ) 

represents the buckled equilibrium condition, the boundary 
conditions of the problem being satisfied by <§) for any values of the 
A. . 

l 
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3. BOUNDARY CONDITIONS 



For the problem of an infinite plate with a hole, two boundaries 
are of concern. At infinity, deflection, slope, and curvature are 
zero. Around the hole, a free edge condition exists. 

The Kirsch solution is used, which already implies zero tractions 
around the hole, and uniaxial tension at infinity, so these conditions 
have already been satisfied. 

For a free edge to exist around the hole means that there is no 
shear, twist, or moment applied at the edge of the hole. Following 
Timoshenko (8), and considering the cylindrical coordinate system 
established in Figure 1, the conditions for a free edge are 



(Q r ) = 0 


(shear = 0) 


(7) 




(twist = 0) 


(8) 


(M r ) =o 


(moment = 0) 


(9) 



Kirchoff (18) has shown that the twist condition can be reduced 
to a statically equivalent shear force. The original requirements 
of zero shear and zero twist can be reduced to one requirement such 
that the sum of the shear forces and the static equivalent of the 
twisting moments equals zero. The conditions then become 



Q y * -p J — O (Kirchoff shear = 0) 

M r = O (moment = 0) 



where Q , 
r 

Qr= -t>£ 

M r = -D 



M^, and M^_ are expressed as 



gV i 3W i 3* W 

Br 1 r ' dir + r x ' a#>- 



a 3 tv v n aw , f d 2 w | 
ac* \r dr + r* a©* j 



M yt = D(i--O) 



r 



a 2 w 

Br^© 



-L . iw. 

r a * d& 



( 10 ) 

(ID 

( 12 ) 

(13) 

(14) 
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When the appropriate differentiations and substitions are made, 
the boundary condition equations become 



s?w i tw i aw (2--0) s*w (3-d) a» a w 

3 r* r dr z r z 'dr + r* '§73©^“ r* 'dQ z 



(Kirchoff 

Shear = 0) 



= O 

y-i 



±ldw,± 

ar 2 - r 1 dr r ‘ a© 2 



= 0 

r=i 



(Moment = 0) 



|W = c£w 
™ 9r 



= o 

Jr=oo 



(15) 

(16) 

(17) 



As previously noted, each member of the family of functions 
must satisfy the boundary conditions in the Ray leigh-Ritz method. 
For a choice of ^ let 



w 



N -1 M =0 L J I 



A nm COS M 9 

0)r 

B nm sin M6> 



(18) 



where A and B„, are arbitrary coefficients that are varied to 
NM NM 

minimize strain energy in the plate. 

The boundary conditions of eq. 17 are satisfied by this series 

since at r = oo an d all of its derivatives are zero. Conditions of 

equations 15 and 16 are satisfied by substituting equation 18 into 

these two boundary condition equations (see Appendix A). Such a 

substitution results in two equations, 15 and 16, with two unknowns, 

^NM ^NM* ^° r ea °k con bination of N and M. A unique value of 

* 

and Y >TW is found which allows each term of equation 18 to 
NM NM n 

satisfy the free edge boundary conditions for any N and M chosen. 

From Figures 1 and 3, it is seen that there is symmetry with 
respect to both the x and y axes. This means that when W is written 
as a Fourier series in 0 (eq. 18), the series contains either cosine 



* A term in this sense means the full expression 
mT V v iFAnmCosMS 

-N j + Xnm ^ Inm • ov 

L y y x J Bmmsiw Me 
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terms or sine terms, but not both. Additionally, the requirement 
of symmetry around both the x and y axes means that only odd or 
even terms in M need be taken. Thus the Fourier portion of 
equation 18 reduces to: 

cos M9, with M = 1,3,5,... 
or 

A x _, cos M0, with M = 0,2,4,... 

NM 

or 

sin M0, with M = 1,3,5, .. . 
or 

sin M9, with M = 2,4,6,... 

This gives four separate expressions for W, one for each 
possible condition of symmetry. The full Fourier series could be 
taken, but if this is done, three quarters of the coefficients, A^ 
or B^. are zero and considerable computer storage space is wasted. 
It is shown later than each of these series represents a different 
buckling mode. 
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4. ENERGY IN THE PLATE 



To this point, a family of functions W(r,0) (equation 18) have 

been developed, which satisfy the boundary conditions independently 

of the coefficients, A._, or B-__, and are adequate to represent 

5 NM NM 

possible shapes of the plate. Attention is now turned upon formu- 
lating an expression for the total strain energy, and then minimizing 
this expression to obtain the critical buckling stress. 

The energy in the plate can be divided into two forms. First, 
there is the membrane energy which is the energy stored because the 
mid-plane is stretched. Second, there is the bending and twisting 
energy which results from the plate buckling out of its plane. 

A basic assumption is that the plate deflection, W(x,y), is 
infinitesimally small, and that the mid-plane strains may be re- 
presented by 



+ 

3| 8 

II 


i law] 


2 \9x 


__ 3r 


i I aw 




2 UtJ 



Y" — 3u, 9r , 3W 3W 

ay 3* d% * ay 



(19) 



Terms of higher order in W than these are neglected, since W 
is infinitesimal; those terms which are retained are essential for 
the calculation of bending energy. 

In the above expressions J— etc., represent the stretching 
of the plate in its own plane. The terms ^ | ^ etc., 

represent the difference between a length of flat plate, and the 
corresponding straight line length when the plate is bent (difference 
between arc and chord) . 
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Stresses in the x and y directions are composed of 
(a) Membrane stresses of finite magnitude, which come 



from the Kirsch Solution and are denoted by (T^ and 
'T X y below; and (b) Infinitesimal bending and twisting 
stresses which do not appear explicitly below, but which 
are involved in the bending and twisting energy. 

The normal forces per unit length (N^,lO can determined from 



N = , etc. The membrane energy in the plate is 

x 

u mem = 2 J/( N x e x + M a + N* 3 rackety 



( 20 ) 



Substituting the strains (eqns. 19), and letting 



II 



MEM 



* U, + U r 



where 




( 21 ) 



( 22 ) 



and 

14 = 




r . 



9W 



dx dtj 



(23) 



an expression for U,„, w is found in terms of the stresses, strains, 

MEM 

and the deflection. 

It can be seen that the membrane energy is composed of two 

terms. One, U^, relates the work done by the stresses, , etc, 

3u 

through a distance , etc. The other term, U 0 , relates the 

O'K l 

I / 3W 

work done by the stresses, 0“ % , etc., through a distance — ( J 5 
etc., as the plate bends out of plane. 

At this point, the assumption is made that the Kirsch 
solution is still valid, even if the plate has buckled infinitesimally. 
This means that is independent of the lateral deflection, W, and can 
be treated as a constant. This assumption is justified by considering 
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that the problem is one of incipient buckling, not postbuckling. 
This assumption considerably simplifies the problem, because it 
uncouples the membrane stress distribution from the buckled 



deflection. 



An expression for the bending and twisting energy, U^, is 
obtained by considering the moments in the plate and the respective 
rotations at a point as functions of plate deflection. Timoshenko (9) 
shows this expression to be; 




aw aw 

dx*- dtf 






9*W _ 3V _ ( 3*W f 
dX x Sy 3 - \3xayl 



dx dy 



(24) 



Adopting the notation Wy-= ) Wg= etc., and transforming 

the previous energy equations to cylindrical coordinates, with the 
proper limits of integration for the problem of a hole in an infinite 
plate, equations 22, 23, and 24 become, respectively 



tL 



U, 



u. 



CONSTANT 



rZTT rOO 2 

j (^w r a + + gL Tto | ^ r de 

*0 J i 




(25) 

(26) 



(27) 

d<9 



CTy-jd© and 'TyQ are obtained from the Kirsch solution (eqns. 1,2, and 3). 
W and its derivatives are obtained from the assumed series (eq 18 ). 
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5. SOLUTION OF THE PROBLEM 



With the expression for total energy (U fc = + \)^ + U^) 

formulated, the next step is to substitute the Kirsch solution, 
and the assumed series for deflection and its derivatives. Then 
this expression must be minimized with respect to each of the 
unknown coefficients, or B . The substitutions and inte- 

grations are shown in Appendices B, C, and D. It is sufficient 
for our purpose here to say that what arises from the integration 
of the expressions for and U^, are two quadruple series, with 
indices I and K summed over the range of N, and J and L summed 

over the range of M, The quadruple series stems from the multi- 

y- 

plication of two double sums. Equations 26 and 27 now become 



U 



2 



U 



3 




I J K L 



(28) 

(29) 



where the flexural rigidity of the plate, D, equals 



Elf 



Q and P are different for each value of I, J, K, or L. 

1JKL 1JKL 

S is the stress applied to the plate at infinity and is factored 
from the Kirsch solution after those expressions are substituted 
into the energy integrals. A and A are the unknown coefficients 
in the series expression for deflection (eq. 18). 



Following the Ray leigh-Ritz procedure, U is differentiated 
with respect to each of the unknown coefficients, so that the 



The unknown coefficients for the cosine series, A^, only will be 

used henceforth, although these could be the sine series coefficients, 

B xtx ., equally as well. 

NM 
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following set of simultaneous equations result ( 



0=^= 


®mkl + Akl 

_K L 1 


+ 


Jk Qnn + £ P m 


An 




11 

£|< 

tola* 

n 

o 


^ ^ (T § PfZKL .) AkL 


+ 


4* *5(2(2 PziZ 


Aia 


> (30) 



} 

The last term in each equation arises from the fact that 

j) = 

As a practical matter, the equations arrived at in this 

investigation are too complicated to permit considering the 

use of a doubly infinite set of terms in the double series. 

Accordingly, it is henceforth presumed that truncated series 

are used and that there is but a finite number of unknown A^'s. 

NM 

Simplification in notation is made at this point if the 
quadruple subscripts are replaced by a double subscript arrange- 
ment. 



Let p> equal the number of terms taken in M, and oC equal 
the number of terns taken in N. Define 

i'=d -op +j 

J'=(K-I)p + L 

Using these equations a unique pair of subscripts can be changed 



to a unique single subscript, so that A^. T = A^ . , , 
etc . 



and Q IJKL Q I'J” 



For example, if 4 terms were taken in M, and we wished to 
convert Q to scme doubly subscripted Q^ lTl , then I = 7, J = 4, 
K = 9 , and L = 2 ; and 



I' = (7-1) 4+4, 



J' = (9-1) 4+2 



= 28 



= 34 



so that Q 7492 becomes Q( 2 8)(34)‘ Similarl y 
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A„ A, 




— > 


Qn 


An. "*■ A 2. 


Qau 




Qai 


Aft -► 


Cl IU2 


-> 


Q u 


etc. 


etc. 







Also, let 

J = (<*-l)-p+p = o<-p 

where J is the largest value of I* and J f respectively. 

Substituting 

pa ..£^— 

^ /2 (I --0 1 ) 

and 

Rf ff * - Qi'J* ar>d Tr'j' = PrV if IVJ' 

RiV = 2*QtV = 2* f^'j* if 

equation 30 simplifies to 




(31) 




9A 5 




Eh 3 - 




J 



When the substitution 

A Eh 2 

is made, and the above set of equations expressed in matrix 



(32) 



notation, the following matrix eigenvalue problem results 

(RX + I)A=0 (33) 

where 
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R,. R 

Ri, R 



12 




R 



(34) 





To *Tia • 

T* r zz . 




T 



(35) 




and 



r 



A, 

Ai 



A 



(36) 



The matrix T is positive definite, since it arose from the 
bending energy integral. The R matrix arose from the membrane 
energy integral, which contains terms from the Kirsch solution. 

This matrix is not necessarily positive definite. Both matrices 
are symmetric. These properties are important to the method 
utilized in solving the matrix eigenvalue problem (see Appendix D ) . 

The question of whether differentiating and setting these 
derivatives equal to zero, produces a minimum or a maximum, has not 
been discussed. Either situation can occur in the problem as it 
has been developed. Some of the eigenvalues obtained in a solution 
of equation 33 represent maxima which are not relevant to the plate 
buckling problem. The criterion for rejecting such answers is the 
set of equations 
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> o 



&u t 

dt\\ 



( 37 ) 



=> o 

BA1 

If the set of equations 37 is satisfied, then is minimized. 

When the set of equations 31 is again differentiated with 
respect to the appropriate coefficient, the set of equations 37 
become 



1 ^= XRi.+ T n -o 
~ X Rzz't' 122 ^ ® 



\ 



(38) 



XR f3 - =» o 

* / 
where an ^ T t , are t ^ le dia § ona l elements of the matrices 

R and T> respectively. All of these terms are positive numbers, 

since they arise from the ’’squares" in the energy integrals (i.e., 

no cross product terms). 

By rearranging equation 38 so that it becanes 



A =■- 3* , 



R 3'3' 



(39) 



It is seen that for some negative values of X , and all positive 
values, the solution of equation 33 represents a minimization of 
energy. 

We are interested in only positive values, so that for any 
answer of interest, a minimization is assured. It should be noted 
that the negative values of A that are admissible stem from 
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buckling modes in the plate that correspond to an applied compressive 



loading. These modes arise from the reversal of sign in the Kirsch 
solution and the associated reversed stress distribution around the 
hole. Modes of this type should be expected from a theoretical 
standpoint . 

When a nontrivial solution has been found for equation 33, 
giving an eigenvalue X and a corresponding eigenvector A, then the 
buckling stress may be found by use of equation 32 and the buckled 
mode shape from equation 18. 
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6. RESULTS AND CONCLUSIONS 



The results of this thesis may be presented compactly in the 
form of a formula for incipient buckling, derived from equation 37, 
viz . 

S cr - Ke(A] (40) 

where S^ is the value of S, the tension at infinity, which causes 
incipient buckling, t is plate thickness, a is the hole radius and 
E is Young's modulus of elasticity. Values of the buckling coef- 
ficient, K, are given in table 1. 

Values of K were obtained from a digital computer solution (see 

Appendix F), for the particular case 
£ 

E = 30 x 10 psi (Young's Modulus) 

-i) = .3 (Poisson's ratio) 

h = .001 hole radii 

Equation 40 permits prediction of for other values of t/a 
and E. However, the effect of different values of -t) cannot be studied 
easily since -0 appears in the bending energy integral and in the 
boundary conditions. All results given are f or i) = .3, but the 
computer program can be modified by changing card listed as 0011 
on page 84 . 

Four distinct mode shapes are apparent, as determined by 
different buckling coefficients appearing for each. A tabulation 
of K is shown in Table 1 below. 
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Table 1 



Values of the Buckling Coefficient, K 



MODE 



K 



1 



7.65 



2 



7.80 



3 



11.7 



4 



12.3 



Since the energy method involves the use of a series which 
cannot be made infinite in the present method, the actual buckling 
coefficients obtained will always be higher than the true value. 
This is because the deflection of the plate is only approximated, 
and the energy never quite minimized. A solution with less plate 
energy could be found, if more terms were taken in the series for 
deflection. 

Different values of the buckling coefficient, K, actually 
calculated by use of the digital computer program are shown in 
Table 2, and are plotted against the number of terms taken in 
both the N and M parameters of the deflection series, equation 18. 
These plots are shown in Figures 4a through 4e. The projected 
values listed in Table 1 are taken graphically from these figures 
as the asymptotes which the curves for buckling coefficient appear 
to approach, when the number of terms in N and M increase. These 
projected values are considered to be the best evaluation of the 
buckling coefficient. 

Figure 4f is taken from the experimental data of Danis (14). 
All of his points for a finite plate are plotted as K versus d/b, 
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Table 2 



Values of the Buckling Coefficient, K, for Various 





Numbers of Terms in the 


Series 






Mode 1 




Terms 


> in M 










4 


5 


6 


7 


8 


9 


6 








11.258 






8 


9.9897 


8.6493 


8.4420 


8.4173 


8.4157 


8.4153 


a 10 


9.1333 


8.0673 


7.8967 


7.8763 


7.8748 


7.8746 


c 12 


8.9990 


7.9743 


7.8060 


7.7857 


7.7840 




14 


8.9927 


7.9683 


7.7987 


7.7777 






g 16 


8.9750 


7.9537 


7.7847 








3 18 


8.9347 


7.9237 










H 20 


8.8897 


7.8913 










Mode 2 




Terms 


> in M 










4 


5 


6 


7 


8 


9 


8 


9.1250 


8.5877 


8.5107 


8.5033 


8.5026 


8.5026 


a 10 


8.5133 


8.0630 


7.9977 


7.9913 


7.9906 


7.9906 


c 12 


8.4310 


7.9860 


7.9197 


7.9130 


7.9127 




‘ H 14 


8.4253 


7.9773 


7.9100 


7.9030 






1 16 


8.4023 


7.9577 


7.8907 








aJ 18 


8.3637 


7.9257 










H 20 


8.3253 


7.8940 










Mode 3 




Terms 


in M 










4 


5 


6 


7 


8 


9 


8 




15.592 


14.545 


14.371 


14.348 


14.347 


a 10 




13.327 


12.616 


12.501 


12.487 


12.486 


c 12 


16.288 


12.679 


12.046 




11.929 




14 


15.959 


12.516 


11.897 


11.793 






S§ 16 


15.910 


12.493 


11.872 








a 18 


15.907 


12.488 










H 20 


15.882 


12.470 










Mode 4 




Terms 


in M 










4 


5 


6 


7 


8 


9 


8 




16.418 


15.765 


15.660 


15.647 


15 . 646 


10 




13.714 


13.318 


13.258 


13.252 


13.252 


a 12 


14.761 


12.929 


12.589 




12.533 




•5 14 


14.463 


12.719 


12.388 


12.337 






w 16 


14.415 


12.682 


12.350 








u 18 


14.412 


12.677 










H 20 


14.392 


12.662 
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the ratio of hole diameter to plate width. The buckling coefficients 
in Table 1, for modes 1 and 2, are plotted on the y axis (for zero d/b) . 

The values for modes 1 and 2 are in fair agreement with experi- 
mental data. No correlation between the theoretical prediction and 
the experimental observation is attempted for modes 3 and 4. 

The buckled deflection of the plate for modes 1 through 4 is 
shown in Figures 5 through 8 respectively. It is seen in these 
figures that the four possible conditions of symmetry are satisfied 
by the four modes. This result is given in Table 3. 



Table 3 

Resultant Symmetry for the Various Mode Shapes 



MODE 



SYMMETRY WITH RESPECT 
TO X AXIS 



SYMMETRY WITH RESPECT 
TO Y AXIS 



1 

2 

3 

4 



f(y) = f(-y) 
f(y) = f(-y) 
f(y) = -f(-y) 
f(y) = -f(-y) 



f(x) = f(-x) 
f(x) = -f(-x) 
f(x) = f(-x) 
f(x) = -f(-x) 



In looking at Figure 4&, it is apparent that mode 2 is 
converging to a value of buckling coefficient that is somewhat 
higher than that for mode 1. This can be qualitatively justified 
by considering the symmetry of buckled shape, and the equivalent 
one dimensional problem. 

Consider two beams, buckled by some unspecified loading, into 
shapes as shown in Figure 9. 




Figure 9 
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It is known from beam theory that the shape of Figure 9b 
has a higher critical buckling load than the shape of Figure 9a. 
This is exactly the relationship between modes 1 and 2 of the 
buckled plate, as can be seen in Figures 5 and 6, respectively. 

There is a small area near the intersection of the hole 
edge and the y axis, which assumes cross sectional shapes similar 
to that in Figure 9a for mode 1, and in Figure 9b for mode 2. 

It is only in these small areas connecting the "top" of the plate 
to the "bottom”, that this sort of bending occurs, so that the 
disparity between the buckling coefficients of the first and 
second plate modes is much smaller than that for the beam modes. 

This same argument applies to the difference between modes 
3 and 4, and to the gross differences between modes 3 and 4, and 
.odes 1 and 2. 
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FIGURE 4a - Plot of the Ducklins Coefficient 
K vs the number of terms taken in 
Mfor mode 1. 
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FIGURE lj.b - Plot of the Buckling Coefficient K 
vs the number of terms taken in M, 
for mode 2. 
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FIGURE 4c - Plot of the Buckling Coefficient K 
vs the number of terms taken in N, 
for modes 1 and 2 
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FIGURE i|d - Plot of the Buckling Coefficient 
K vs the number of terms taken in 
N for mode 3. 
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Plot of the Buckling Coefficient K 
vs the number of terms taken in N, 
for mode 1+. 
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FIGURE 4f 
PLOT OF K vs d/b 
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FIG 5 - MAP OF DEFLECTION FOR MODE I 
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FIG 8 - MAP OF DEFLECTION FOR MODE 4 
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APPENDIX A 



SATISFACTION OF THE BOUNDARY CONDITIONS 
The free edge boundary conditions for the hole are stated 
in equations 15 and 16 as 
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(A- 1) 



(A- 2) 



The plate deflection, W, is taken as 
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where the coefficients A__. and B ATN . are varied to minimize the plate 

NM NM 

energy, and the coefficients and Y are chosen so that the 

boundary conditions expressed by equations A-l and A-2 are satisfied. 

The derivatives of W are tabulated in Table 4. When these 
derivatives are substituted into equations A-l and A-2, and evaluated 
at r = 1, the following equa tionsresult 
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where 



S, = N(N + 0 - N-0- M 2 -0 (A- 6) 

S a = (N-H)fN + 2)-(N+l)-0- M 2 t) (A-7 ) 

S 3 = (N+2)(N+3)-(N+2)t)-M 2 -0 (A- 8) 

t , = - M(N+l)(N+2) -I- N(N*l) + N + (2- -0)M 2 N 4- (3-t>)m 2 (A- 9) 

t 2 =-(N+ , )(N+2)[N+^)+(N+l)(N+2)+-(N+>) + (2--0)M 2 fNfl)+ O-OJM 2- (A- 10) 

t 3 =-(N+2)(N4-3)(Wf4) + (N+2)(N43) HW+2f)+ (Z--d) M^N+2) + (3-t))M* (A-ll) 



Equations A-4 and A-5 must hold for all values of 9, which can 
only be true if 



S, 4* S z X NM 4- S 5 — O 

ti + tz^NtA + ^ ^ 



(A- 12) 
(A- 13) 



for any value of N and M. Thus, a value of and Y^. can be 
found for each value of N and M, such that the boundary conditions 
are satisfied. 

By application of Cramer's Rule: 
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The above results are incorporated into SUBROUTINE BOUND (see 
APPENDIX F) which generates values of and for given values 
of N and M. 
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APPENDIX B 



EVALUATION OF THE ENERGY INTEGRALS 



The total energy in the plate, from equations 25, 26, and 27, is: 

U t — U| +■ U 2 + U3 

where 
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U 2 = ~ / [cTy-W* + 2^0 + r dn d9 



U 3 = l 
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dr d0 



Note that these integrals contain second degree terms in the 
derivatives of the deflection, W. Since W is represented by a 
series, care must be taken to insure that the series are properly 
multiplied. 

The product of two double sums, both summed over the same range 
for the subscripts N and M, is a quadruple series, with dummy subscripts 
I and K summed over the range of N, and K and L, summed over the range 
of M. 
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It is seen that the cross product terms above do not combine 

directly into one term. One cross product term involves r ^ the other 

-K 

involves r . For this reason, cross product terms in the energy 
integrals B-l and B-2 must be identified and separated. 

As an aid in keeping the cross product terms separated, the 
notation is adopted that starred (*) functions are functions of I 
and J, and that non-starred functions are functions of K and L, in 
place of N and M, respectively. Equation B-l may be written 




+ 4 < Ww^W* + W r W e *)]rdrd 



(B-3) 



and with the substitution of the Kirsch solution becomes 
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Equation B-2 is now 
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where 
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A listing of the expansion of the products of the derivatives 
necessary for the evaluation of Ug and are shown in Tables 5 
and 6, respectively • Further presentation of the energy integrals is 
in the notation delineated in those tables „ 

Equations B-4 and B-5 now become 
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When the integration, in € is done, equation^ B-9 and B-10 reduce to 
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Tables 5 and 6 



Tabulation of the Integrals of the Radius Dependent Terms 
NOTATION : 

G i> G 2 >...E^, E , etc. are functions of I, J, K, and L. 

x means X^ as used previously, only as a function of K & L. 

x* means X xtx , as used previously, only as a function of I and J. 

NM 

y means as used previously, only as a function of K and L. 

y means Y^ as used previously, only as a function of I and J. 

THE INTEGRALS ARE TABULATED ON THE FOLLOWING PAGES. THE FIRST 
NINE TERMS IN EACH COLUMN ARE SUMMED TO FORM THE WHOLE EXPRESSION. 
THE TERMS a(J,L) or b(J,L) SHOULD BE CHOSEN FOR THE UNKNOWN 
COEFFICIENTS or B ffl RESPECTIVELY. 
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Table 5 



Tabulation of the Compacted Farm of the Terms in the 
Membrane Energy Integral 
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Table 5> (continued) 
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Table 5 (continued) 
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a(J,L) = cos J© cos L© 
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Table 5 (continued) 
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Table 5 (continued) 
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JLx 


4 . 


JLx 


'(I+K+3 ) 


T 


TT+K+5T 


JLxx" 


4. 


JLxx" 


(i+k+1 ).T 


T 


( I +K+ 6 ) 




4. 


JLxy* 


(l+K+5) 


T 


(I+K+7) 


JLj 


4. 


JLy 


(I+K+I 4 J 




( I+K+ 6 ) 




4. 


JLyx* 


(l+K+5^ 


T 


(I+K+7T 


JLyy* . 


4 . 


JJyy" 


( i+k+ 6 ) 


T 


(Y+K+5T 



where 

= sin J9 sin L9 a(J,L) = sin J9 sin L9 

= cos J9 cos L© b(J,L) = cos J© cos L© 
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Table £ (continued) 



f 



7a 



JK 

(I+K) 




LI 

Ii+kT 



J(K+l)x 

- rrmfr 



L(I+1 )x“ 
(i+K+iy 



J(K+2) 
“ (I+K+2 






L( 1+2 )y“ 

■■(I+f+2T 



JKx*' : " 

( I+K+1 ) 

J(K+1 )xx" 
Tl+K+2) 



J(K+2)yx 
' (I+K+3) 

JKy" 

" Tt'+K+2')'' 

J(K+1 )xy“ 
~ (I+K+3) " 



- W 



where 

a(J,L) = -sin J9 cos LG 
b(J, L) = sin LG cos JO 



LIx 

(I+K+1) 

L( 1+1 )xx“ 
'TI+K+2T* 



LIy 

- (i+ifer 

L(I+1 )yx* 

“ (I+K+3) 

. l(i+ 2 )yy — 

(I+K+ljj 

where 

a(J,L) = -sin LG cos JO 
b(J,L) = sin JO cos LG 



L( 1+2 ) 



U+2jxy 

rnwfr 
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Table 5 (continued) 



* 

s 

1 

II 

c 


„ LI 


oa \ i+k+2 ) 


ob " " TT+^+2) 


J(K+l)x 

u+^+3')' 


L( 1+1 )x“ 
(l+K+3) 


J(K+2 )y 

(I+K+I 4.5 


L(I+2)y*" 

•n+k+ii.T 


JKx" 


LIx 


TT+K+3 ) 


(l+K+3) 


J(K+i )xx" 
"Tl+K-Hj. ) 


L( 1+1 )xx* 
(i+K+ij.) 


J(K+2)yx* 

‘ T 


L(I+2)xv* 

Ti+k+5) 


JKy* 


LIy 


TI+k-Hj.) 


Ti+k+Ij. ) 


J(K+l)xy“ 

Ti+k+5) 


L( 1+1 )yx* 

" "(T-lS?)' ' 


J ( K+2 ) yy*" 
"(l+K+6) 


L( 1+2 )yy* 
' ("I’h-K+6 ) 


where 


where 


a(J,L) = -sin J© cos L© 


a(J,L) = - sin : 


b(J, L) = sin L© cos J© 


b ( J,L) = sin J© 
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Table 5 (continued) 




where 
a(J,L) 
b( J,L) 



JK 

" ( I+K+I4. ) 



J ( K+l ) x 
( I+K+5 ) 

J(K+2 )y 
(l+K+6) 



JK x* 

( I+K+5 ) 

J(K+1 )xx“ 
TT+k+SI — 



J(K+2)yx 
( I+K+7) 



(I+K+6) 



J(K+1 )xy‘ 
'(I+K+7y 

J (K+2)yy' 

TT^TBT" 



= -sin JO cos L0 
= sin L9 cos J9 




where 
a( J,L) 
b(J,L) 



LI 

" (I+K+IO 
L(I+l)x* 

(I+k+5) 




LIx 

(T+KTcry 

L( 1+1 )xx“ 

'(Y+K+5T“ 

L ( I +2 ) xy* 

-umr 




L(I+1 )yx -> 
(I+K+7) 



= -sin L0 cos J9 
= sin J0 cos L9 
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Table 6 



Tabulation of the Compacted Form of the Terms in 
the Bending Energy Integral 



-oo 



G, = 



W W" 
rr rr 



dr 




‘ A U , ' A KL 

b( J,L) 




co, 

W’w 



r rr 



- o 



dr 



} a ( J,L) • A t t *A. 



= S 



2a 



or 



IJ KL 



|b( J,L) 'B I j'B kl 




oo 

W W* 
r rr 



dr 



s 2b 



"a(J»L)-A 

or 

_b{ J,L) 




f W rr»69 



dr 




a ( J, L ) ' A j j * | 

or 

_b( J,L) 'Bjj'Bj^ J 




\ W W* 

I rr 99 

i * 2 



dr 




a( J,L)' 
or 

_b( J,L) 'Bjj'B^j 



G, = 



oo 

r w w* 

j r r 



I 



dr 






a( J,L)> A IJ 'A KL 'i 



or 



_b( J,L) • b ij‘ b klJ 






OO 

w ee w r 



J- 



dr 



,a( J,L) 'Ajj’AjCL. 




oo 

W 



w” 
©Or 



dr 



" S 5b 



| a( J,L) • Aj j 'Aj^l 
or 

_b( J,L) ' b ij ,b kl 
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Table 6 (continued) 



OO 



°6 ■ 



w ©e w ee 



7T 



dr 



g 6 * 



a(J,L) 

or 

b ( J,L) 



oo 



W W“ 
r© rO 



dr 



= g 



7 ' 



a(J,L). 

or 

b(J,L). 



G 8a = 



OO 



' Vre 



dr 



} a ( J , L ) • 
1 or 
jb(J,L) 



s 8a* L,°. r 



G 8b = 



oo 



! W-W" 

/ Q r© 



J r ' 



- dr 



j a. ( L ) 



= 



co 



V« 



1 



dr 



= 6c 



a( J,L)* 
or 

b( J,L)* 



* a ij' a kl 

* b ij b klJ 

• a ij' a kl' 

•ht-hasl 

J 

a ij a kl1 
b ij b klJ 

a ij* a kl" 
b ij V. 

a ij ,a kl*I 

B I J ’ B Kxj 
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Table 6 (continued) 



S 1 



IK( 1+1 ) (K+l) 
(I+K+21 



( 1+1 ) (K+l ) ( 1+2 ) (K+2 )xx* 

+ 

+ (I^)(K+2Hl+^)(K+3jj7^ 



K(K+l)(I+l)(I+2)x* 

(I+K+3") 



1(1+1) (K+l) (K+2 )x 

+ (TTETT) 

. K(K+ 1 ) (1+21(1+3)/' 

(I+K+I 4 .) 

+ KI+l) (K+2 j )(K+3)y 

. ( 1+1 ) ( 1+2 ) (K+ 2 ) (K+3 )yx* 

(I+K+5) 

. (K+l ) (K+2 ) ( 1+2 ) ( 1+3 )xy* 

(I+K+5 ) 



where 

a(J,L) = cos J9 cos LG 
b(J,L) = sin J9 sin L9 
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Table 6 (continued) 



_ _ IK (K+l ) 

g 2a " (l+K+2) 


IK( 1+1 ) 
g 2b ” (I+K+2) 


KK+l) (K+2)x 

" "TitK+31 


K(I+1) (I+2)x* 

- (i+k+5T“ — 


_ I(K+2) (K+3)y 

(I+K+4) 


K( 1+2) (1+3)/' 

- u+fr+Cr 


(1+1 ) (K+l )Kx* 
" (I+K+3) 


(K+l )I(I+1 )x 
(I+K+3) 


(1+1) (K+l ) (K+2 )xx* 

■ — cr+K+in 


(K+l) (1+1) (1+2) xx* 

- — ci+K+rn 


(I+l)(K+2)(K+3)yx* 

- (T+K+5T "" 


(K+l) (1+2) (1+3 )xy* 
(I+K+4) 


- 


(K+2)I(I+l)y 

(I+K+4) 


(1+2) (K+l ) (K+2 )xy" 
( I+K+5 ) 


(K+2 ) ( 1+1 ) ( 1+2 )yx* 

(i+k+5) 


( 1+2 ) (K+2 ) (K+3 )yy* 
( I+K+6 ) 


(K+2) (1+2) (1+3 )yy" 

(i+K+sy 


where 


where 


a(J,L) = °os J0 cos L0 


a(J,L) = cos J0 cos L0 


b(J,L) = sin J© sin L© 


b(J,L) = sin J0 sin L0 
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Table 6 (continued) 



S 3a = 



I ( 1+1 ]Lr 

(I+K+2) 



KM)ri 

(I+k+3) 



I ( I+l ) L' 



i+l JL y 



(I+l) (I+2)L 2 x* 

(l+K+3) 

( I+K+4. ) 
(I+l)(I+2)L 2 yx* 

(T+ftt^r — 



(I+l) (I+2)L 2 xx* 



(I+2)(I+3)L 2 t* 

(ti-R %T"- 



g 



3b 



K ( K+l ) J 2 
* (I+K+2) 

K(K+1) J 2 ^ 

nm+TT 



# 



2* 



K ( K+l ) J^y 
(K+1)(K+2)J 2 x 

-'irmyr — 

( K+l ) ( K+2 ) J 2 xx" 

— — 



(K+2) (K+3) J^yx*' 
(I+K+5) 



where 

a(J,L) = cos J© cos L© 
b(J,L) = sin J© sin L© 



where 

a(J,L) = cos J© cos L© 
b(J, L) = sin J© sin L© 
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Table 6 (continued) 



II 

t5* 


IK 


g 5a ~ 


J^L 


TT+K+z) 


(I+K+27 


+ 


K(I+l)x* 

TTttTT 


+ 


J^(K+1 )x 

wmr ■ 


+ 


K(I+2 )y* 
(I+K+llI 


+ 


J 2 (K+2)y 

(l+K4fr 


+ 


( K+l ) lx 

(Y'+Y+fT 


+ 


J^Kx* 

(I+K+3T 


4- 


(K+l ) (1+1 )xx" 


4 . 


J 2 (K+l)xx* 




(I+K+ 1 +) 


T 


(I+K+I 4 .) 


+ 


(K+l ) (I+2)xy* 

Ti+k+^) 


+ 


J 2 (K+2)yx"’ 

(i+k+$T 


+ 


(K+2)Iy 

TT+t+rrr 


+ 


Ay* 

(I+A) 


+ 


(K+2 ) (I+l)yx* 

— iT+K+^r — 


+ 


^(K+Dxy* 

rr+K+n 


+ 


<K+2)(I+2Wj* 

(l+K+ 6 ) 


+ 


J 2 (K+2)yy* 

(i+k +£1 


where 




where 




a(J,L) 


= cos J© cos L© 


a(J,L) 


= cos J© cos 


b( J,L) 


= sin J© sin L0 


b( J,L) 


= sin J© sin 
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Table 6 (continued) 



= l 2 i 

g 5b li+i+2) 

+ l 2 (i+i}** 
"l i+k+3) 

+ 4W 1 

. L 2 Ix 
(I+K+3) 

. L 2 (I+1)xx' : 

(I+K+4) 



nr+K+5 ) 
, l 2 i 



r+fan 



, L 2 (I+1 )yx 

+ 



. L 2 (I+2)yy ; 

+ (i+fi+SI^ 1 



g 6 



J 2 ! 2 

(I+K+2) 

J 2 ! 2 ** 

(i+k+3) 



5T 



jW 

(I+K+3) 

jVxx* 
( i+K+i|.) 




J 2 ! 2 






TT+K+ 

fogy 

(I+K+ljJ 

J 2 L 2 yx* 
(i+K+f)'' 

fo 2 yy * 

( I +K+o J 



where where 

a(J,L) = cos J© cos L9 a(J,L) = cos J© cos L© 

b(J,L) = sin J© sin L© b(J,L) = sin J© sin L© 
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Table 6 (continued) 



_ IJKL 
s 7 (l+K+2 ) 



g 8a 



JK 

(I+K+2) 



IJ(K+l)Lx 
(I+K+3 ) 



J(K+l)x 

IW 






(I+K+lj.) 



J(K+2)y 
( I +K+I4. ) 



, (1+1) JKLx* 

+ TT+in3T“ 



JKx~ :: 

(I+K+3)' 



(1+1) J(K+l)Lxx“ 
( I+K+i (. ) 



J(K+1 )xx*' : 

nw“ 



(I+l) J(K+2)Lyx* 
(I+K+5) 



J(Ki2)yx 

(I+K+5) 






(1+2) JKLy* 
(I+K+k) 



JKy* 

(I+K-fi).) 



(1+2) J(K+1 )Lxy* 

+ rm^r 



J(K+l)xy” 

(i+3+fr 



( 1+2 ) J( K+2 ) Lyy^ 
(I+K+M 




where 

a(J,L) = sin J© sin L© 
b(J,L) = cos J© cos L0 



where 

a(J,L) = sin J© sin L© 
b(J,L) = cos J© cos L© 
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Table 6 (continued) 



LI 


g 9 


JL 


K 8a TT+K+2) 


( i+k!+2 ) 


L( 1+1 )x“ 

- ■ TY + K + '3r 


+ 


JLx* 

TT+k+TT 


(I+K+Ij.) 


+ 


JLy* 

TT+k^) 


LIx 

( I+K+3 ) 


+ 


JLx 

T7T&5) 


L(I+1 )xx* 

- "rr+Tfjr 


+ 


JLxx* 

(I+K+ij.) 


L( 1+2 )xt* 

- (i^tr 


+ 


TT * 

JL xj 
(I+K+£) 


LIy 

( I+X+k ) 


+ 


JLy 


L( I+l )yx" 

- TT^rr” 


+ 


__ * 

JLyx 

'(i+K+fTF 


L(I+2)yy“ 

“ 1 Ti+wT 


+ 


JLyy* 

(I+K+&) 


where 


where 




a(J,L) = sin J9 sin L© 


a(J,L) 


= sin 


b(J,L) = cos J9 cos L© 


b( J,L) 


= cos 
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APPENDIX C 



INTEGRATION OF SINE-COSINE TRIPLE PRODUCTS 



When evaluating the energy integrals for the plate, it is 
necessary to evaluate the following two integrals 
-27T 



cos A9 cos B@ cos C© d© 



and 



( 



r Z7T 

I s ~ I (sin AO sin B 9 cos CB) de 



(C-1) 



(C-2) 



where A, B, and C are integers. 

The above integrals can best be evaluated by expanding the 
integrands using the identities 



^[siN(*x+y)+ SiN(Tt-y)] 


(C-3) 


l[siN(x+y)-SiN('x.-y)] 


(C-4) 


rj^cos (-x+y)*cos(x-y)] 


(C-5) 


r|^cos(7C-y)- cos(a+y)j 


(C-6) 



The object is to reduce the integrals to the form 



I 



27 r 

sin D9 dQ 



or 



1 



27 r 

COS DQ d& 



(C-7) 



0 -'o 

where D is an integer sum or difference of A, B, and C. Note that the 
sine integral in C-7 is always zero and that the cosine integral is 
nonzero only if D is zero. 

Expanding equation C-2, using identity C-3, we get 



r znf r 

J $ = 1 Sin A©[sin(B 4-C)©+- SiNfB-c)© 



dQ 



(C-8) 
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Continuing the process with identity C-6 



I $ = 4 j [cos(A-e-C)6-cos(A+B+c)&4.cos(A-b+c)0-cos(A<-B-c)©jd© 



I g is now expressed as the sum of four integrals 



Is ~ I| ^ I2 ^I^ + 14 

where 



I 



1 



^/2 if (A-B-C) = 0 
0 otherwise 



-'11/2 if (A-ffi-fC) = 0 
0 otherwise 

'IT/2 if (A-B-fC) = 0 
0 otherwise 
- 1*/ 2 if (A+B-C) = 0 
0 otherwise 



A similar expansion for I c (equation C-l) can be made 
results are the same as above, except that the signs of I 
are positive. 



(C-9) 
(C- 10) 



The 

and I. 
4 



76 . 



APPENDIX D 



SOLUTION OF THE EIGENVALUE PROBLEM (Rx + T) A = 0 



The canonical form of the matrix eigenvalue problem is 



(ix-b)x = o 



(D-l) 



where I is the identity matrix, X an unknown eigenvector, \ is 
an unknown scalar parameter (eigenvalue) , and B a known coefficient 
matrix. 

Equation D-l has a solution if determinant of 



If this determinant is expanded, a polynomial equation in h results. 
However, if the determinant is of order greater than three or four, 
it is generally inefficient to obtain this polynomial equation and 
then solve it. 

There are alternative procedures for solution of the eigenvalue 
problem. A method which gives extremely accurate eigenvectors and 
eigenvalues for symmetric matrices is one named for Jacobi, which 
is presented by Fox (1). Brown and White (13) give a comparison of 
the accuracy of the Jacobi method with that of other methods available. 
The necessity for extremely accurate eigenvectors is shown later. 

Consider the eigenvalue problem (equation 33) 



where both R and T are symmetric matrices and T is positive definite. 
Dividing equation D-3 by X , and making the substitution 




(D-2) 




(D-3) 




(D-4) 



gives 




(D-5) 
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Any symmetric matrix with real elements can be reduced to a 
product of the eigenvector and eigenvalue matrices, so that 

I = SRa T (d-6) 

where g. is a diagonal matrix of the eigenvalues (d^) of T, and Q 
is a matrix of normalized eigenvectors (q^) , obtained from a 
solution of the problem 

(i<lt-l)cu=0 (D-7 ) 

Note that ^ T = q} since the matrix T is symmetric and has orthogonal 
eigenvectors . 

Equation D-5 may now be written 

(HE£ T y “B)a = ° (d-8) 

Premultiplying by , we get 

(g§ T| j“$ T E)A = ° (D-9) 

T 

Factoring a Q out of A = I A = (£ ^ A, and letting 

V = Q r A (D- 10) 

the problem reduces to 

(ey-| T i^)Y = ° (d- ii) 

Since is positive definite, D contains all positive terms. 

* % ” % 

By again premul tiplying^ this time by D , factoring D , and making 
the substitution 

2 = d' / s -Y (d- 12) 

equation D-ll reduces to 

(xy“g Va ‘i T lS£' /i )i =0 (D-i3> 
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and becomes 

(l«3-§}2 = o (d- 14) 

which is of canonical form. 

% -% 

Note that D and D are easily obtained by taking the square 
roots, and the square root reciprocals, respectively, of the diagonal 
terms. Also note that B is symmetric. The problem is now well 
suited for solution by use of the Jacobi method. It is seen that 
inaccuracies in obtaining eigenvectors would cause compounded 
errors in the final solution for A and \ . 

The problem is returned to original form with the substitutions 

X = - ^ (D- 15) 

and 

A = Qp' /a Z (D- 16) 
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APPENDIX E 



SUGGESTIONS FOR EXTENSION OF THE PROBLEM 

There are two obvious directions for extension of this problem 
of a thin plate with a circular hole. The buckling problem for 
a strip of finite width and the postbuckling problem for both finite 
and infinite plates remain to be treated analytically. 

The postbuckling problem is the more difficult of the two 
because of the coupling between the stress distribution and the 
buckled deflection. Johns (2) outlines the postbuckling situation 
for plates and suggests several methods of solution. Solution of 
the postbuckling problem would be important because of the revised 
stress concentration factors that would be obtained; cf. Houghton 
and Rothwell (16). 

Extension of this problem to a finite strip, of width 2b, would 
require the following: 

a. Using the Howland solution (17) instead of that due to Kirsch, 
to obtain the stress distribution in the plate. 

b. Satisfying additional boundary conditions for two "free" 
straight edges. 

Integration of the energy terms for the finite strip is between 

k 

the limits r = d and r = — — - , so that the integration with 

Sin q & 

respect to 0 becomes more complicated than for the infinite plate. 

This integration is over the interval 0 to 27T. Wylie (10) 
presents a method for the evaluation of real trigonometric integrals 
over this interval. The method involves substitution of appropriate 
complex variables for the sine and cosine terms, and then evaluation 
of the residues inside a unit circle. 
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There are two possible methods for satisfaction of the 
additional boundary conditions of the finite strip. Introduction 
of additional unknowns, similar to and in equation A- 3, 

into the assumed series for deflection could be made. These 
unknowns could be determined by simultaneous solution of all the 
boundary condition equations. Another method, the Lagrangian 
multiplier method, introduces dummy variables so that additional 
restraining equations, those expressing the boundary conditions, 
can be solved simultaneously with the energy expressions. 
Budiansky and Hu (12) present a clear example of the use of this 
method in minimizing the plate energy, while satisfying plate 
boundary conditions . 
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APPENDIX F 



COMPUTER PROGRAM LISTINGS 

The following pages contain a listing of the computer programs 
used in the solution of the problem. All programs are written in 
FORTRAN IV for use on an IBM System 360 computer. 

The control card listing on the first page is effective as of 
September 1967. Users of this program should check current procedures. 

The maximum terms to be taken in N and M, and the plate thickness 
•read into the computer by data card (see comment cards regarding 
input requirements). Young* s modulus, E, and Poisson's ratio are 
established by cards serial 0010 and 0011 in the main program. 

Printed output from the computer is a column matrix EIVU which 
is an array of all of the values of S in psi for the parameters 
of plate thickness, Poisson's ratio, Young's modulus, terms in N, 
and terms in M. Matrix EIVU is followed by a map of the buckled 
plate deflection for the minimum value of tensile critical buckling 
stress , 
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